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ABSTRACT 



Context. The generation of X-rays in accreting T Tauri stars (TTS) is thought to be predominantly due to energy dissipation in mag- 
netic fields, but alternative X-ray generation mechanisms have been proposed, such as heating in accretion shocks near the stellar 
surface, or in shocks forming in jets. 

Aims. We test models and trends discussed in the literature using X-ray data from the classical TTS T Tau. 

Methods. High-resolution spectroscopy from the Reflection Grating Spectrometers on XMM-Newton is used to infer electron densi- 
ties, element abundances and the thermal structure of the X-ray source. We also discuss the ultraviolet light curve obtained by the 
Optical Monitor, and complementary ground-based photometry. A high-resolution image from Chandra constrains contributions from 
the two companions of T Tau N. 

Results. The X-ray grating spectrum is rich in emission lines, but shows an unusual mixture of features from very hot (as 30 MK) 
and very cool (1-3 MK) plasma, both emitted by similar amounts of emission measure. The cool plasma confirms the picture of a soft 
excess in the form of an enhanced O vn/O vm Lya- flux ratio, similar to that previously reported for other accreting TTS. Diagnostics 
from lines formed by this plasma indicate low electron densities (5: 10 10 cm -3 ). The Ne/Fe abundance ratio is consistent with a trend 
in pre-main sequence stars in which this ratio depends on spectral type, but not on accretion. 

Conclusions. On the basis of line density diagnostics, we conclude that the density of the cool "soft-excess" plasma is orders of 
magnitude below that predicted for an accretion shock, assuming previously determined accretion rates of (3 - 6) x 10~ 8 M o yr _1 . 
We argue that loading of magnetic field lines with infalling material suppresses the heating process in a part of the corona. We thus 
suggest that the X-ray production of T Tau is influenced by the accretion process although the X-rays may not form in the bulk of the 
accretion footpoints. 

Key words. Stars: coronae - Stars: formation - Stars: pre-main sequence - X-rays: stars - Stars: individual: T Tau 



1. Introduction 

Low-mass pre-main sequence stars maintain magnetic fields 
that are thought to provide the energy for high-energy par- 
ticles and radiation as shown by the ir strong X-ray and ra- 
dio e mission (for rev i ews, s e e, e.g. , Feigelso n & Montmerld 
1999; FeigelsonetaD 120061: iGudell 120021 l2004l) .~ Vigorous 
sources of X-ray radiation have been found among clas s 
I protostars (e.g., [Kovama et al.l [19941 llmanishi et al.l 200ll). 



accre t ing classical T Tau stars (C TTS, e.g., Walter & Kuhi 
ll98ltlFeig elson & DeCamrjli 198 JJ), and weak-line T Tau stars 
(WTTS; e.g.. IWalteiHl986t iNeuhiiuseret al.lll995l) that are no 
longer surrounded by a thick, accreting disk. A number of ob- 
servational features have supported a picture in which X-ray 
emission is due to magnetic coronal processes similar to mech- 
anisms seen in the solar corona: rotational modulation points 
to a c oncentration of hot plas ma in various magnetic active re- 
gions dFlaccomio et al.L 12005); flares heat plasma on time scales 
of minutes to hours, reaching temperatures of a few tens of 
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MK dWolk et all |2005[): and electron temperatures continuousl y 
exceed 10 MK (e.g.. llmanishi et alj|2001t ISkinner et al]|2003l) : 
such temperatures cannot be achieved by liberating gravitational 
energy from infall, nor by shocks in moving gas (e.g., jets) com- 
monly observed around young stars at other wavelengths. Also, 
high-energy electrons observed through their gyrosynchrotron 
emission give direct evidence for particle acceleration, for which 
magnetic reconnection processes provide an ideal basis. 

Nevertheless, accreting pre-main sequence stars add circum- 
stellar disks and accretion flows to the picture of the stellar envi- 
ronment; stellar magnetic fields may reach out to the inner bor- 
der of the disk from where they guide li ghtly ionized mass to the 
stellar surface dCalvet & GullbringHl998l) . WTTS and CTTS are 
therefore ideal objects to test to what extent the stellar environ- 
ment controls coronal magnetic fields and energy release. Except 
for disks and mass accretion streams, the two types of stars ap- 
pear to be identical in nature, occupying the same region in the 
Hertzsprung-Russell Diagram (HRD). 

Evidence for some sort of influence of the disk environ- 
ment onto X-ray production is broad: i) Star-disk magnetic fields 
may brake stellar rotation, thus damping the magnetic dynamo 
and inducing weaker overall magnetic activity. In the Taurus 
Molecular Cloud, for example, CTTS rotate on average less 
rapidly than WTTS (characteristic rotation periods of P a 8 d 
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and 2—4 d, respectively; Bou vier et al.l 1993). Many T Tau stars 
may, however, emit at the empirical sat uration limit, Lx/^boi ~ 
10 -3 ' 5 , for rotation periods below 10 d dPreibisch et all 120051) . 
so that a difference between the dynamo operation of CTTS 
and WTTS may be difficult to identify, ii) Despite possible sat- 
uration effects, it has now been clearly established that both 
the X-ray luminosity, Lx, and Lx/^boi are, for otherwise sim- 
ilar stellar chara cteristics, lower by a factor of « 2 for CTTS 
than for WTTS dNeuhauser et all 1 19951: iPreibisch et all l2005t 
iTelleschi et all l2006al) . It is unclear in what sense accretion 
is responsible for this result, iii) On the other hand, the elec- 
tron temperature of the X-ray emitting plasma may be higher 
in CTTS than in WT TS dlmanishi et alll200ltlTsuiimoto et all 
120021 ITelleschi et all |2006a). iv) Specific changes in the time 
behavior of X-rays were also observed in outburst phenomena 
ascribed to strong accretion events in FU Ori-type and EX Lupi- 
type variables (e.g., a strong flux increase of hard emission dur- 
ing outburst, iKastner et alj 12004a, 120061; o r a softening of the 
plasma during outburst, lAudard et al.ll2005l) . 

Recently, a more direct role of accretion has been sug- 
gested for the generation of X-rays in CTTS, namely emis- 
sion from gas that has been shock-heated in the accretion 
flow near the stellar surface. While free-fall velocities of a 
few 100 km s _I on CTTS are sufficient to heat shocks to a 
few MK, a problem is X-ray absorption by the gas stream 
itself, and by the surrou nding photosphere i f the shock is 
buried below the s urface dCalvet & Gullbrinj, 1 19981: iLamzinl 
1 999: iDrakel 12005). However, the temperatures attained in such 
shocks are clearly not sufficient to explain the extremely hot 
X-ray sources of CTTS with temperatures exceeding 10 MK. 
On the other hand, in the exceptional case of the CTTS TW 
Hya, X-ray spectroscopic eviden ce for a predominantly cool 
plasma led IKastner et al] (2002) to suggest that X-ray pro- 
duction by accretion may be relevant in this star. This con- 
jecture has been further supported by evidence for unusu- 
ally high electron densities of n e > 10 11 cur 3 derived from 
density-sensitive line flux ratios of Ovn and Neix in high- 
resolution X-ra y spectra of t he CTTS TW Hya, BP Tau , 
and V4046 Sgr dKastner et all 120021: IStelzer & Schmittl r2 004; 
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2006), supported by flux ratio s of Fexvn lines dNess et al 



2005). IStelzer& Schmittl d2004l) further argued that high Ne/Fe 



and N/Fe abundance ratios may point to accreting gas that is 
depleted of Fe, presumably because Fe (and several other ele- 
ments) easily condense onto dust grains further out in the disk, 
leaving an accreting gas phase enriched by Ne and N. 

Two large survey projects of X-ray emission in star- 
forming regions hav e addressed the above questions anew. 
Preibisch et al.l d2005l) found no support for an accretion sce- 
nario in in the Chandra Orion Ultradeep Project (COUP); the 
bulk of the observed X-ray emission is related to very hot, coro- 
nal plasma. Accretion may in fact suppress the X-ray luminos- 
ity by mass loading of magnetic fields. High-resolution X-ray 
spectroscopy collected in the framework of the XMM-Newton 
Extended Survey of the Tauru s Molecular Cloud (XE ST) has 
shown a series of new facets dTelleschi et all l2006bl) : i) very 
high (k, 10 11 cm 3 ) electron densities are not recorded in all ac- 
creting pre-main sequence stars, posing problems for the accre- 
tion scenario even in the most favorable case of 100% accretion 
surface filling, ii) Abundance "anomalies", in particular a high 
Ne/Fe abundance ratio, are common to many T Tau stars, regard- 
less of whether or not they are accreting. On the other hand, the 
Ne/Fe abundance ratio seems to depend on the the stellar spec- 
tral type, iii) While the bulk plasma observed in CTTS shows 



very similar thermal properties as plasma in WTTS, there is soft 
excess emission in accreting stars; this excess is visible as (and is 
defined by) an enhanced ratio of the O vn/O vra Lya line fluxes 
compared with WTTS or zero-age main-sequence (ZAMS) stars. 
Because the maximum formation temperatures of O vn and O vm 
are 2 MK and 3 MK, respectively (and O vm is formed up to con- 
siderably higher temperatures), an enhanced ratio indicates the 
presence of additional plasma below 3 MK. Observation of this 
excess requires high-resolution spectroscopy that resolves emis- 
sion lines formed at 1-4 MK such as lines from O vn or N vi and 
N vii. 

To further our understanding of these issues, detailed case 
studies are needed. Ideally, a bright source will produce detailed 
high-resolution grating spectra of the line-dominated range be- 
low 1 keV. This allows us to assess anomalies in the thermal 
structure and electron densities from line ratios. To access lines 
formed at low (1^ MK) temperatures (located predominantly 
at wavelengths longward of 17 A), the photoelectric absorption 
should be moderate (hydrogen column density A^h of no more 
than a few times 10 21 cm -2 ). Only few CTTS fulfill these re- 
quirements. One of the most outstanding among them is the 
CTTS T Tau N. We report detailed X-ray spectroscopic and 
imaging studies of this object, also searching for emission re- 
lated to its two embedded companions, T Tau Sab. 

The specific goal of this paper is to use T Tau as a test case 
for hypotheses of X-ray production related to magnetic and ac- 
cretion processes. We add complementary information from op- 
tical and ultraviolet observations. The plan of the paper is as 
follows. In Sect. |2j we introduce the T Tau system, and describe 
our observations in Sect. [3] Sect.|4]presents our imaging studies, 
Sect.|5]introduces light curves, and Sect.|6]presents results from 
our spectroscopic analysis. Sect.|7]contains a discussion and our 
conclusions. 

2. Target 

Historically, the optically revealed northern component, T Tau 
N, was the prototype of accreting low-mass pre-main sequence 
stars. Subsequently, a strongly variable IR companion was found 
ft/7 to its south, now designated as T Tau S. It was soon 
recognized t o be a protostar probably co-eval with T Tau N 
dDvck et all [l982) and dominating the system flux at A > 
2/im. The southern component was later foun d to be d ouble , 
with a separa t ion of only ft' 09 or 13 AU dKoreskol 120001; 
iDuchene et all 120021) . An orbital analysis of astrometric ra- 
dio and infrared data infers masses of (2.73 ± 0.31) M Q and 
(0.61 +0.17) M Q for Tau Sa and Sb, respectively dDuchene et all 
2006), suggesting that the former is an embedded Herbig Ae star. 
The mass of T Tau N is also high er than the masses of typical 
T Tau stars, namely M * 2 4 M n dGiidel et all l2006al) . and its 
rotation period is P = 2.8 d dHerbst et al. , 1986). The distance to 
T Tau is * 140 pc dLoinard et alll2005h . 

The T Tau triple may be at a critical transition stage between 
protostars and fully optically revealed T Tau stars; some of its 
observational properties appear to be related to its particular ori- 
entation that places us into an opt imal position for observations , 
van Lang evelde et alJdl994al) and lvan Langevelde et alJ dT994b) 
found evidence for infall of molecular mate rial and a col l imated 
outflow, pointing at 13° to the observer. ISolf & Bohml (1999) 
suggest a pole-on view for T Tau N, which is also clearly sup- 
ported by mm observat ions of a well-deve loped, near-face-on 
massive accretion disk dAkeson et alj. 1 19980 . The latter authors 
point out that the disk is similar in mass to the solar nebula in 
which our planetary system formed. The pole-on view of the 
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Table 1. X-ray observing log and coordinates of T Tau 



Parameter 


XMM-Newton 


Chandra 


XEST number 


01-045 




ObsID 


0301500101 


3732 


Start time (UT) 


2005-08-15 13:52:13 


2003-10-24 06:11:52 


End time (UT) 


2005-08-16 12:55:22 


2003-10-24 13:37:02 


Exposure time 


82989 s 


26710 s 


Coordinates:" 






Boresight RA 


04h 21m 59.4s 


04h 21m 59.69s 


Boresight decl. 


19deg 32' 06" 


19deg 32' 2U! 6 


X-ray: RA 


04h 21m 59.44s 


04h 21m 59.42s 


X-ray: decl. 


19deg 32' 05'.' 8 


19deg 32' 06! 6 


T Tau N: RA* 


04h 21m 59.43s 


04h 21m 59.43s 


T Tau N: decl.* 


19deg 32' 06'.' 4 


19deg 32' 06'.'4 



a All coordinates are J2000.0 

b For origin of coordinates, see Smit h et all d2003h 



N component results in a relatively low extinction, while from 
larger angles all T Tau components would appear to be em- 
bedded accreting young stars. A CO emission ring, identified 
with the outer walls of a bi-conical outflow region, provides ev- 
idence that we are looking down into a low-extinction cavity 
which the stellar win d has nearly swept free of molecular gas 
dMomose et al.L[T996h . T Tau Sa a ppears to be sur r ounded by a n 
opaque disk seen nearly edge-on dDuchene et all I2005L [2006). 
The disk then is no more than 3 AU in radius given the close 
orbit of the Sb companion (with a semi-major axis of 11.6 AU 
and a periastron distance of * 6 AU). 

The T Tau system has been detected as a radio source with 
properties of an extended wind but also of non-therm al radio 
emission (Phil lips et aIlll993l:[Skinner & Browrill994l) . includ- 
ing gyrosynchrotron emission from an outflow (Ray et al. 1997) 
and strongly polarized emission from magnetic fields inferred t o 
be of kGauss strength dSmith et all 120031: lLoinard et all 120051) . 
The T Tau system was well detected in X-rays wit h the ROSAT 
PSPC (Neuhaus eretalL fT99l IStdzeret al.. 2000), with a qui- 
escent count rate of (0.03+0.009) cts s _1 , and multiple energetic 
and rapid flares developing on ~1 hour timescales and reaching 
up to 0.2 cts s . 



3. Observations and data reduction 

Table Q] summarizes basic parameters of our XMM-Newton and 
Chandra X-ray observations, including the centroid coordinates 
of the X-ray images of T Tau. Both exposures comprise data 
of four stellar/substellar systems that are thought to be mem- 
bers of the Taurus Molecular Cloud: these are, apart from the 
T Tau system, the M0 CTTS IRAS 04187+1927, the WTTS 
RX J0422.1 + 1934, and the brown dwarf 2M J04221332+1934. 
All have been detected as X-ray so urces. Furthe r details on the 
former two stars are summariz ed in Giidel et al. (2006a), while 
the brown dwarf is described in lGrosso et all (120061) . 

We obtained an *83 ks exposure of t he T Tau sys t em on 
August 15-16, 2005, with XMM-Newton dJansen et all 1200 ll) 
using standard observing setupsQ Specifically, the European 
Photon Imaging Cameras (EPIC) obtained images with an angu- 
lar resolution of a 4 - 5" (FWHM), using the full windo w with 
a diameter of 30' for the PN camera (Striide r et all 1200 lb . while 
the MOS1 and MOS2 cameras dTurner et all 12001) observed 



in the large-window (W3) and small-window (W2) configura- 
tions that constrain the central CCD to a field of 5. 5 'x5. 5 'and 
1.8'xl.8', respectively. These constrained window modes were 
chosen to prevent the detectors from pile-up effects in case of 
strong flares (but no strong flares that would induce pile-up were 
actually observed). All EPIC cameras used the medium blocking 
filter. 

The two Reflec tion Grating Spectrometers (RGS; 
Iden Herder et all [20011) were operating in spectroscopy mode, 
providing X-ray spectra with a resolution of approximately 
60-70 mA (FWHM) in the wavelength range of 5-35 A. Two 
chip failures earlier in the mission left a gap in the RGS1 
spectrum in the range of 10.6-13.8 A and in the RGS2 spectrum 
in the range of 20.0-24. 1 A. 

The Optical Monitor (OM; iMason et al.ll200lT) observed in 
imaging and high-time resolution mode with the UVW1 ultra- 
violet filter inserted. The latter transmits light roughly in the 
2500-3500 A range. The imaging field of view is a square with 
a size of 17'xl7'. The high-time resolution recording applies, 
however, only to a central window of 10.5"xl0.5" in which T 
Tau was positioned. 

The X-ray data ana lysis of the XMM-N ewton observations 
is described in detail by Giidel et al. (2006a). For details on the 
OM data reduction, see lAudard et al. (2006). The data reduction 
made use of standard SAS v.6.1.0 software (i.e., the epchain, 
emchain, rgsproc, omfchain, and omichain tasks). Spectra and 
images were extracted with the evselect task. For the EPIC data, 
we used canned response matrices appropriate for SAS v.6.1.0, 
whereas ancillary response files (ARFs) were generated using 
the SAS task arfgen. The RGS response matrix was generated 
with the rgsrmfgen task. Background radiation fluctuated be- 
tween a high and a low level during the observation, each oc- 
cupying approximately half of the time. The high particle radia- 
tion can affect accurate background subtraction in both RGS, in 
particular in long-wavelength regions where the observed, true 
continuum is low and suppression of the source spectrum by 
photoelectric absorption may be severe. We have experimented 
with either extracting only low-background data, or with con- 
straining the width of the source extraction region in the RGS 
detector plane. We found that extracting 70% (down from stan- 
dard 90%) of the PSF width was optimum, and this choice was 
preferred to selecting in time while keeping the extraction width 
at the standard value. 

The nearby WTTS RX J0422. 1 + 1934 was bright enough 
(amounting to 41% of the cou nt rate of T Tau in the EPIC de- 
tectors, see iGiidel etai1l2006al) that it slightly contaminates the 
background region on the RGS detector. The source region of 
this star was therefore excluded from the background region of 
T Tau. It does not directly contaminate the T Tau spectrum, how- 
ever, as it is located at the edge of the RGS detector, far outside 
the standard spectral extraction width for T Tau. 

We further obtai ned a Chandra X-Ray Observatory 
dWeisskopf et all 120021) exposure of 26 ks duration, u sing the 
High Resolution Camera (HRC-I, iMurray et all l2000h with a 
spatial resolution of m 0.5" but no energy resolution in the X- 
ray range. The HRC field of view is about 30' in diameter. The 
data were reduced using standard CIAO v.3.0.2a software toolfl 
New level 1 events were created with the hrc_process_events 
task, and level 2 events after applying filters and good time in- 
terval (GTI) selection. A light curve analysis revealed no signif- 
icant variability. A formal source detection was performed with 



1 see XMM-Newton User's Handbook (http://xmm.vilspa.esa.es/ 1 



2 http://cxc.harvard.edu/ciao/guides/hrc_data.html 
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the wavdetect algorithm in order to obtain precise centroid coor- 
dinates for the T Tau source. 

We obtained ground-based optical observations of T Tau dur- 
ing three nights around the Chandra observations. Photometric 
monitoring was performed in the UBVR bands with the 48-cm 
AZT-14 reflector at the Mt. Maidanak Observatory during the 
nights of 2003 October 22/23, 23/24, and 24/25. For these ob- 
servations we used a single-channel photoelectric photometer 
in pulse-counting mode, and the photometric data were reduced 
to Johnson's international system. We observed HD 27560 as a 
comparison star and HD 284415 as check star. Other ob serva- 
tional information can be found in iMel' niko v et al. (200]|). 



4. A high-resolution X-ray image of T Tau 

The T Tau triple with a separation of 077 between the N and 
Sab components cannot be resolved by the XMM-Newton EPIC 
cameras. The Chandra HRC image, on the other hand, could 
marginally reveal binary structure if the N and S components 
showed comparable flux levels. The centroid coordinate of the 
stellar X-ray image (Table[TJ is very close to T Tau N. 

Fig. Q] shows the HRC-I image of the T Tau source. We 
modeled a point-spread function (PSF) using the CIAO task 
mkpsf. For this, we used a circle with a radius of 15 pixels 
(with a pixel size of 0.13") around the centroid. The PSF was 
then fitted to T Tau image in the Chandra sherpa software 
dFreeman et all 12001 1) (the version implemented in CIAO ver- 
sion 3.3.0.1 was used). The model also contained a spatially 
constant background. We found the centroid to be slightly but 
insignificantly offset from the true stellar position (offset "X-ray 
minus true" of -0711 in RA and 0725 in declination^). 

Next, two PSFs were fitted, now assuming two sources with 
fixed (offset-corrected) positions, the correct relative coordi- 
nates, and the same PSF width. Although a very faint extension 
toward the S is suggested in the Chandra X-ray image, the best- 
fit amplitude at this position is only (6.5+2.7)% of the amplitude 
of the T Tau N model, and is therefore at best a marginal detec- 
tion (at the 2.4cr level). 

In summary, the Chandra HRC image (Fig. [TJ shows that 
T Tau N is the principal contributor to the X-ray source. This 
finding is supported by spectral analysis below that shows low 
photoelectric absorption compatible only with T Tau N. A 
small contribution from T Tau S cannot, however, be excluded. 
Because its visual extinction is Ay > 8 mag (for T Tau Sb, 
iDuchene et al.ll2002h . we expect a hydrogen column density of 
Nr > 1.6 x 10 22 c m" 2 , assuming standard gas-to-dust ratios 
dVuong et all 120031) . Some counts could therefore be expected 
at energies above « 1.5-2 keV only. We estimated an X- 
ray luminosity of T Tau S with the PIMMS softwar^ adopt- 
ing A^h = 1.6 x 10 22 crrT 2 , an electron temperature of either 
« 7 MK or 28 MK (corresponding to either of the higher tem- 
peratures found for T Tau N), and 6.1% of the total counts. 
We found L x ,o.3-io keV = 4 x 10 30 erg s" 1 and L x ,o.3-io kev = 
1.4 x 10 30 erg s _1 , respectively, for the two temperatures. The 
true values may be higher if Nh is higher, but lower if we over- 
estimated the count rate contributions in the marginal fit. In any 
case, given the small contribution compared to T Tau N and lack 
of strong significance, we will assume that the observed XMM- 
Newton spectra are from T Tau N. 



3 The 90% source location error circle in Chandra has a radius of 
about 0.5", see Chandra Proposers' Observatory Guide v.8. 

4 http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 
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Fig. 1. The T Tau source on the Chandra HRC-I detector (pixel 
size: 07 1 3). The positions of T Tau N and Sab are marked by the 
crosses. The white circle around the expected T Tau N position 
indicates the 90% source location error for Chandra. There is 
a slight systematic offset of about 0727 of the X-ray image rel- 
ative to the true position of T Tau N (-0711 in RA, +0725 in 
declination in the sense "X-ray minus true"), which has not been 
corrected for in this figure. 



5. The optical, UV, and X-ray light curves 

Fig. [2] shows the X-ray and ultraviolet light curves from the 
XMM-Newton observations. 

We extracted the combined EPIC PN and MOS1 
background-subtracted light curves for the softer and harder por- 
tion of the spectrum (energies of 0.4-1.0 keV and 1.0-7.3 keV, 
respectively; MOS2 was not used because of difficulties in 
defining a background curve in the small-window mode). These 
light curves are shown in the upper two panels of Fig. [2] Both 
count rates are slowly decaying, probably from a flare with a 
peak that appears to have been recorded at the beginning of the 
observation. The harder flux decays more rapidly but then also 
levels off earlier, leading to a hardness increase in the second 
half of the observation, illustrated in the third panel that shows 
the ratio "hard:soft". The fourth panel shows the light curve for 
the entire spectral range. Note that during the first w 1350 s, 
the PN camera did not record data (interval before the time 
marked by the left arrow). We scaled the MOS1 light curve 
by an average factor to the estimated PN+MOS1 count rate. 
Similarly, only the PN observed during the last 2830 s (time 
after the second arrow). The light curve has again been scaled 
from the PN count rate to the expected total count rate. 

The OM UVW1 light curve is shown in the bottom panel. It 
was recorded in a series of exposures, each with a duration of 
~ 3540 s, followed by a gap of 340 s. A few exposures were not 
successfully recorded. The high count rate induced large coinci- 
dence losses which led to inaccurate calibration of the incident 
count rate in some exposures, visible as obvious discontinuities 
with respect to the count rates in previous or following expo- 
sures. The affected exposures were manually corrected by mul- 
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Fig. 2. Light curves of the T Tau source, extracted from the 
XMM-Newton observations. From top to bottom: Combined 
EPIC PN+MOS1 light curve for soft (0.4-1 keV) counts; hard 
(1-7.3 keV) counts; hardness ratio, i.e., ratio between hard and 
soft count rates; total light curve (0.4-7.3 keV; the arrows mark 
the times between which both EPIC cameras were operational); 
and the OM UVW1 light curve (dots illustrate bin count rates 
before manual correction; see text for details). Horizontal bars 
give time intervals for each OM exposure. 



tiplying their count rates to produce a smoothly connected light 
curve (Fig. |2 the uncorrected fluxes are also shown; the count 
rates in the first exposure were not corrected and are uncertain 
given the long time gap before the second exposure). 

Three features are to be mentioned: i) There is no detailed 
correlation between the X-ray and the UV light curves, ii) Both 
curves show, however, a general decay in the first half of the 
observation, perhaps as a consequence of the weak flare seen in 
X-rays. The time of the UV peak is uncertain, given the unknown 
correction factor for the first exposure. It may have occurred be- 
fore the start of the observations, iii) The UV light curve in- 
creases in the second half of the observation, in contrast to the 
X-rays. However, the X-ray hardness also increases, and both 
UV flux and hardness show a peak around t — 6 x 10 4 s. It 
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Fig. 3. Optical UBVR light curves of T Tau, collected during 
three consecutive nights (2003 October 22-25). The abscissae 
give JD-2452900.0, and the ordinates are in magnitudes. Both 
day-to-day variations and short-term fluctuations (in the U band) 
can be seen. 



is possible that both features relate to a flare-like heating event 
that produced both hot coronal plasma (seen in the harder X-ray 
band) and enhanced emission in the chromosphere and transition 
region (seen in the UVW1 band). 

The light curves in UBVR obtained at Mt. Maidanak 
Observatory during the ground-based monitoring period in 2003 
around the observing campaign with Chandra are presented in 
Fig. [3] No overlap with the Chandra observing interval was 
achieved. These optical obse rvations near l y cove r the 2.8-day 
photometric period of T Tau ( He rbst et al.L I1986I) which corre- 
sponds to the rotational period of the star. Our data show day -by- 
day variability of T Tau: while on October 22/23 the brightness 
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Fig. 4. Statistical analysis of variability in the U band light curve 
of T Tau (top) and the check star HD 284415 (bottom). The solid 
histograms give the number of measurements with a given mag- 
nitude offset from the average (after removing a slow trend for T 
Tau). The dashed lines mark the maximum rms error for a mea- 
surement at the average magnitude. 



of the star is rather stable, the brightness must have risen in all 
bands between the first and the second observation; prominent 
slopes can be seen in the light curves during the next two days, in 
opposite directions, first decreasing in brightness, then increas- 
ing again. Given that the rotation period of the star is 2.8 d, at 
least two minima and two maxima must have occurred during 
the stellar rotation covered by these observations. Magnetically 
active stars indeed often sho w two longitude intervals with pre- 
ferred magnetic activity (e.g., IBerdvugina et al.l ll998). The am- 
plitude of this variability seems to increase with decreasing 
wavelength and therefore reaches its maximum for the U band. 
This is particularly evident during the October 24/25 observa- 
tion when the U-band curve approaches the B-band brightness. 
This day-by-day variability is similar to the long-term variations 
in the XMM-Newton UV light curve. 

Another type of variability which is suggested in the U-band 
curve is defined by fast small-amplitude variations. In order to 
estimate the accuracy of our observations, we calculated his- 
tograms of the brightness distribution for both the check star and 
T Tau itself (Fig. |4j». For the second and third T Tau observation, 
we first removed a slow trend in the brightness, using a polyno- 
mial fit that describes the nightly trend but no shorter features. 
The x-axis gives the deviation of points from the mean magni- 
tude after removal of the trend. The maximum rms error of a 



single measurement in the U band calculated from the observa- 
tions is 0.01 1 mag for the check star and 0.021 mag for T Tau (in 
the U band the check star was roughly three times brighter than 
T Tau). We plot these errors with dashed lines in the histograms. 
The histogram of T Tau suggests that some of the variability is 
due to real fluctuations beyond statistical errors, while this is not 
the case for the check star. 



6. High-resolution X-ray spectroscopy 

We present the combined and fluxed RGS 1 +2 spectrum in Fig. [5] 
Given the significant background radiation level, we prefer to 
first show the full spectrum before background subtraction and 
the background spectrum itself separately. The potential influ- 
ence of background features on the spectrum at the source po- 
sition can thus be examined in detail. The unsubtracted spec- 
trum also retains maximum signal-to-noise ratios to reveal weak 
line features. (We also inspected the spectrum extracted for low- 
background episodes only. Apart from a lower S/N, it very 
closely resembles the spectrum shown here.) For the subse- 
quent quantitative spectral analysis, we will use background- 
subtracted RGS spectra and provide separate illustrations (see 
below). We mark all spectral line features that contained ^ 
(3 - 4) x 10~ 6 photons ctrT 2 s , identifying the dominant ions 
contributing to the flux (a few of these lines may not be individ- 
ually recognized in the spectrum). Dotted vertical lines point to 
the location of a few important fainter lines that remained unde- 
tected. 

The spectrum is rich in lines, but also reveals a strong con- 
tinuum shortward of 17 A. The most unusual aspect of the T 
Tau X-ray spectrum is its mixture of features associated with ex- 
tremely hot plasma, such as the continuum and several lines of 
highly-ionized Fe formed above 10 MK, and lines attributed to 
very cool plasma such as those from N vn and O vn. The ratio 
between the fluxes in the O vn He-like line triplet at 21.6-22.1 A 
and the O vm Lya line at 18.97 A is near unity. This is unusual 
for magnetically active main-sequence stars and WTTS, while 
other CTTS have shown similar ratios, as noted in a survey of T 
Tau X-ray grating spectra presented by Tellesc hT et alj (|2006b). 
Among main-sequence s olar analogs, only ina ctive, old stars re- 
veal such low flux ratios (Tell eschi et al.Ll2005l) . We also note the 
unusually clear appearance of two He/3 lines of Ovn, at 17.7 A 
and at 18.6 A and the presence of the Lya line of N vn at 24.7 A. 
On the other hand, several features in the unsubtracted spectrum 
coincident with the Hear lines of N vi at 28.5-29.5 A and the He/3 
lines of N vi at 23.3 A and at 23.8 A are not significant because 
the background spectrum contains similar features at the respec- 
tive wavelengths. We emphasize that the identified features of 
O vn do not suffer from irregular background. 

We have performed spectral fits in XSPEC dArnaudl 1 1996b . 
using the vapec collisional ionization equilibrium thermal 
line+continuum model combined with a common photoelectric 
absorption component, parameterized by the neutral hydrogen 
column density Nu- We use background-subtracted spectra for 
this analysis. The source spectrum is sufficiently rich in spectral 
features that a determination of several element abundances is 
warranted. We note that alternative methods such as reconstruc- 
tion of emission measure distributions from extracted line fluxes 
could in principle be applied; however, most of the important 
line features required for such an analysis are heavily blended, 
and from no element (including Fe) do we see lines from more 
than two ionization stages at sufficient signal-to-noise ratios. We 
therefore prefer to use template spectra in XSPEC that take into 
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Fig. 5. Fluxed, combined RGS 1 +2 spectrum. The background has not been subtracted, but its spectrum is overplotted as a red, dotted 
histogram, and also separately shown in the narrow, lower panels on the same flux scale. The bin width is 58 mA. The brightest line 
features contributing to the model fit are labeled. A few important features that have not been detected are marked with dotted lines. 



account all unresolved blends and faint lines present in the code, 
while the principal thermal components are determined by the 
strongest lines and the continuum. A comparison between these 
two approaches for similar but brig hter X-ray spectra from solar 
analogs has shown little difference (Tellesc hi et all 12005). 

We simultaneously fitted spectra of RGS1 and RGS2 be- 
tween 7 A and 25 A, complemented with a portion of the MOS 1 
spectrum in the range of 1.5-9.35 A. We use only one MOS 
spectrum confined to the harder range because we prefer to give 
full weight to the high-resolution RGS spectra wherever possi- 



ble. The MOS spectrum complements information on the Mg, 
Si, S, and Fe abundances and constrains the hottest plasma com- 
ponents. All spectra were rebinned so that they contained at least 
15 counts per bin before background subtraction. Three thermal 
components were required for a good fit, defined by their temper- 
atures 2\2,3 and the associated emission measures EMi ; 2,3- The 
results are shown in Table [2] together with 68% and 90% confi- 
dence ranges. The reduced X 1 is excellent (1.09 for 382 de grees 
of freedom = dof). We also list the "average temperature" T av , 
which has been defined as the logarithmic average of the three 
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Table 2. Results from the 3T model fit 



Parameter 




best-fit 


68% range 


907c range 


N H [10 ZI cnT 2 ] 




0.49 


0.46-0.55 


0.43-0.59 


Ti [MK] 




1.76 


1.30-2.01 


0.94-2.16 


T 2 [MK] 




7.42 


7.11-7.69 


6.91-7.87 


T 3 [MK] 




28.8 


27.4-30.6 


26.5-32.3 


EM] [10 53 cm" 3 ] 




9.3 


5.5-17.8 


3.8-25.8 


EM 2 [10 53 cm" 3 ] 




4.4 


3.4-6.2 


2.8-7.7 


EM 3 [10 53 cm" 3 ] 




5.0 


4.6-5.4 


4.1-5.6 


Abundances:" 










N 




0.31 


0.00-2.92 


0.00-5.62 







0.43 


0.28-0.75 


0.28-1.24 


Nc 




0.82 


0.61-1.02 


0.50-1.19 


Ma 






U. JH — U.UJ> 




Si 




0.27 


0.17-0.37 


0.12-0.44 


c 




U. jU 


U.l /— U.44 


n nn n za 


Fe 




0.28 


0.21-0.35 


0.18-0.40 


T m [MK] 




5.2 


4.2-7.3 




)i c [10 1U cm" 3 ] 






< 2.6 


< 8 


^,0.1-10 ^ 




2.02 






L X,0.3-10 [ 1Q31 ^ 


M c 


1.52 






dof 




1.09 

382 







" Rela tive to solar photospheric ab undances of Ander s~& Grevessel 
( 1989) except for Fe, which refers to lGrevesse & Sauvall ( I1999T) 
b for the 0. 1-10 keV range, using a distance of 140 pc 
c for the 0.3-10 keV range, using a distance of 140 pc 



temperatures, weighted with the respective emission measures; 
its 68% range is an approximation calculated using Gaussian 
error propagation. This temperature is considerably lower than 
T av = 12.8 - 17 6 MK f ound in the XEST survey from EPIC 
data (Gudel_etalJ,|2006a). This is because the EPIC spectral fit 
does not recognize the softest emission for which the cool com- 
ponent at 1.8 MK is responsible. If we calculate r av for the two 
hotter components, we obtain T av = 15.7 MK, in agreement with 
the EPIC survey value. The fits are illustrated in Fig. [6] for two 
spectral ranges of the RGS detectors. The RGS data are shown 
at a higher resolution (fewer counts per bin) than actually used 
for the fit, in order to better illustrate the individual line features. 

The unusual mixture of very hot (*30 MK) and very cool 
(«1.8 MK) plasma is confirmed. In the intrinsic, unabsorbed 
model spectrum, the O vn resonance line is in fact the strongest 
line in the entire X-ray spectrum, exceeding even the O vm Lya 
line by a factor of 2.0. Such ratios are otherwise typical only 
for extremely inactiv e stars with cool coronae, such as Procyon 
dRaassen etalll2002l) . 

The amount of cool emission measure is similar to the to- 
tal emission measure of the hotter two components although the 
precise ratio strongly depends on A^; the latter sensitively affects 
the soft part of the spectrum to which the cool plasma mostly 
contributes. The value of Nn ~ 4.9^9g x 1 21 cm" 2 is somewha t 
higher than found from EPIC data alone (Giid el et all |2006a), 
N n = (3.4 + 0.1) x 10 21 cirr 2 (from a 2-T analysis). Part of this 
discrepancy could be induced by slight, systematic calibration 
inaccuracies at the long-wavelength end of the RGS. Further, for 
spectra of limited S/N, an acceptable spectral fit could overesti- 
mate the cool EM if A^h is slightly too high. However, the basic 
result is evident in the flux ratios of the Ovn and Ovm lines 
at similar wavelengths. We find a ratio between the observed 
fluxes in O vn and O vm of 1 .06 + 0.29. From the model, we in- 
fer that about half of the O vm Lya line flux is due to the two 
hotter components. The other half must be due to a cooler com- 
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Fig. 6. Spectral fit to the background-subtracted RGS 1 (black) 
and RGS2 (red) data in the 11-18 A (top) and 18-23 A (bot- 
tom) ranges. Note the missing data in the RGS 1 spectrum below 
13.8 A and in the RGS2 spectrum above 20 A. These gaps are 
due to chip failures earlier in the mission. The lower panels in 
each plot show the residuals in terms of sigmas. 



ponent and therefore amounts to about 50% of the total Ovn 
flux at 21.6-22.1 A. Keeping the hotter two components fixed 
in the model but enforcing A^h = 3.4 x 10 21 cm 2 as in the 
XEST survey (iGiidel et alj|2006al noting that A^ of this mag- 
nitude principally influences temperature and emission measure 
of the coolest component), we calculated conditions for T\ and 
EMi to produce the observed flux ratio, given the contributions 
from the hotter components. There is a family of solutions for 
(7*1, EMi) below T\ =2.1 MK, with the smallest emission mea- 
sure still being 17% larger than EM2, and this solution is ob- 
tained for T\ - 1 .59 MK. We conclude that even for the lower 
A^h, the coolest emission measure is comparable to the interme- 
diate and hotter EMs. 

The element abund ances (given with respect to solar pho- 
tospheric abundances of Anders & Grevesse 1989 except for Fe 
for which we used the value of lGrevesse & S auval 1999 as a ba- 
sis) are common among T Tau stars dTelleschi et al.Ll2006bl) . We 
note a relatively large abundance ratio for Ne/Fe (= 2.9 + 1.0). 
The abundance of N is very poorly determined, given the very 
low S/N of any potentially detected line. The usual trend to- 
ward an inverse first ionization effec t in ma gnetically active stars 
dBrinkman et all 1200 ll IGiidel et all |2001|) is marginally recov- 
ered. 
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We used the Ovn He-like line triplet at 21.6-22.1 A 
to infer electron densities in the coolest plasma compo- 
nent. In high-density environments, the forbidden transition 
(f, ls2s 3 5'i — > Is 2 'So) at 22.1 A can become suppressed 
by electron collisions in favor of the intercombination transi- 
tions ( i, ls2p 3 P L2 -> \s 2 1 S ) at 21.8 A dGabriel & Jordanl 
Il969t iBlumenthal etall 1 19721 see iTelleschi et alJ l2006bl for 
applications to high-resolution spectra of other pre-main se- 
quence stars). The density-sensitive range approximately cov- 
ers 10 ll) cirT 3 n e S; 10 12 cm -3 . The limited S/N makes indi- 
vidual extraction of the three line fluxes difficult. The i line is 
not significantly detected, while the r and / lines are present (at 
the 2.3 and 3.6cr level, respectively), suggesting low densities. 
Formally, the unabsorbed fluxes at the distance of the Earth are 
(5.14 ± 1.41,0.91 ±0.95,3.70+ 1.58) x 1(T 4 photons cm" 2 s" 1 in 
the (r, i, f) lines, resulting in a best-fit // i ratio of 4.0. However, 
we note that there is an additional constraint on the relative line 
fluxes. The so-called G ratio, (r + i)/f, where r, i, and / de- 
note the fluxes in the resonance, intercombination, and forbidden 
lines, is determined by the thermal structure of the plasma. We 
therefore fixed our best-fit thermal model and made use of the 
density calculations implemented in the vmekal code in XSPEC 
to infer the densities. We only fitted « e and, for optimum adjust- 
ment of the line fluxes, the emission measure of the coolest com- 
ponent that predominantly forms the lines. The high f/i ratio 
requires low densities. The best-fit values are below the density- 
sensitive range (in agreement with the above f/i estimate) so 
that we quote only the 68% and 90% upper limits in Table [2] 
These are approximately 2.6 x 10 10 crrT 3 (68% confidence) and 
w 8 x 10 10 cirT 3 (90% confidence) using Gehrels weights on our 
spectral data. 

The X-ray luminosities, Lx, in Tablewar e given for the 0.3- 
10 keV range (adopted for the XEST survey. iGQdel et alj|2006ah 
and also for the 0.1-10 keV range to include important con- 
tributions from the coolest component. T he latter ind e ed dou - 
bles Lx compared to the EPIC analysis in lGiidel et alJ (|2006a). 
Not considering the contributions from the cool plasma, we find 
Lx = 1.0 x 10 31 erg s" 1 for the 0.3-10 keV range, nearly identi- 
cal to the XEST survey value of 0.94 x 10 31 erg s" 1 based on a 
two-component spectral fit to the EPIC PN data. 



7. Discussion and conclusions 

7.1. Observational signatures 

We now discuss the X-ray properties of T Tau N as found in 
this paper in the context of X-ray features that have previously 
been associated with accreting pre-main sequence stars, and with 
properties seen in the U-band and optical light curves. We ad- 
dress, in turn, the following characteristics: i) Flux ratios in He- 
like triplets, indicative of electron densities; ii) cool plasma; iii) 
abundance anomalies; and i v) optic al an d ultraviolet va r iability . 
We refer to ITelleschi et ail d2006bl) and ITelleschi et all d2006c) 
for further detailed discussions of high-resolution X-ray spectra 
from pre-main sequence stars. 



5 We confirmed these estimates using Monte Carlo simulations based 
on the observed source spectrum, the observed background spectrum, 
the same response matrix and exposure time, and the thermal struc- 
ture and abundances as found in our best fit to the entrire spectrum; 
we obtained an 68% limit at 2 x 10 10 cm" 3 and a 90% limit around 
ss 9 x 10 10 cm" 3 . Conversely, if the true density were 10" cm" 3 , we 
found that only one out of ten statistical realizations converges to an 
inferred density below 



He-like line flux ratios and electron densities. High elec- 
tron densities in cool plasma could be a consequence of accre- 
tion shocks. The post-shock temperature in the accretion column 
near the stellar surface is T — 3v 2 fim p /l6k, where the veloc- 
ity v is about 0.9 times the free-fall velocity Vg = (2GM/R) 1 ^ 2 
(ICalvet & GuUbringi Il998h . m p is the proton mass, k is the 
Boltzmann constant, and the mean molecular weight fi « 0.62 
for a fully ionized gas. For T Tau, we use M = 2AM e and R = 
3.6R Q as listed in lGiidel et all (l2006al) . Then, v ff = 500 km s" 1 
and T ~ 2.9 x 10 6 K. The soft excess is thus compatible with 
expected post-shock temperatures of the free-falling gas. 

The accretion rate of T Tau is M acc » (3 - 6) x 10" 8 M o yr" 1 
(IWhite & Ghezll200UlCai"vet et all 120041) . For the accretion lu- 
minosity (neg lecting viscous dissipat ion), we use the expres- 
sion given bv ITelleschi et all d2006bl) : L acc30 ~ 600MM_ 8 /,R, 
where L acc , 30 = L ac c/(10 30 erg s" 1 ), M = M/M e , M acc _ 8 = 
M acc /10~ 8 M o yr" 1 and R = R/R Q . We thus obtain L acc ,3 = 
1200 - 2400, i.e., at least two orders of magnitude more than 
required to produce the radiative output in the softest X-rays. 
Some energy will be conducted to cooler la yers, and the m ajor- 
ity of the X-rays may be entirely absorbed dLamzinj|1999l) . 

Lastly, we can estimate the shock density, using the strong- 
shock condition «2=4«i, where n\ and n 2 are the pre-shock and 
post-shock densities, the former relating to the mass accretion 
rate and the accreting surface fraction F (filling factor - a s sumin; 
that all flows produce similar shocks). ITelleschi et al.l d2006 
gave the expression 



«2 



4x10" M_ 8 

£3/2^1/2 ~ 



[cm -3 ]. 



(1) 



Typic al values for F are F = 0.1 - 10% (ICalvet & Gullbrinjd, 
1998). For T Tau, the density should then be n 2 = (3.6 - 7.2) x 
10 9 M_ 8 /Fcnr 3 =(l.l-2.2)xlO n AFcm- 3 > 1.1 x 10 12 cnT 3 , 
where the latter value is for the extreme case of F = 10%, 
while we obtain n f * (1.1 - 2.2) x 10 14 cm" 3 for F = 0.1%. 
ICalvet et al.l (120041) derived a filling factor of F = 3.1% specif- 
ically for T Tau. This implies n e as (3.5 - 7.1) x 10 12 cm" 3 . 
These values are orders of magnitude larger than the densities 
(or upper limits thereof) inferred from the O vn triplet. In con- 
trast, the T Tau density upper limit (« 8 X 10 10 cm" 3 at the 90% 
confidence level) is fully compatible with the measured densities 
of (mostly magnetic ally active) main-sequence coronae reported 
bv lNess et al. (2004). These are distributed between < 10 10 cm" 3 
and about 10 11 cm" 3 (derived from O vn), the median and mean 
logn e values being 10.36 and 10.37 + 0.06, respectively (as 
derived from the Kaplan-Meier estimator analysis in ASURV 
which also takes upper li mits into account; Fei gelson & Nelson! 
119851; lLa Valley et alj|1992l) . 

We also inspected the Ne ix density-sensitive triplet formed 
at somewhat higher temperatures («4 MK). The signal-to-noise 
ratio is low and the Ne ix triplet is strongly blended with several 
lines of highly ionized Fe. Nevertheless, there is a clear indica- 
tion of a strong / line in that triplet as well, again pointing to 
low densities (the Ne triplet is density-sensitive in the range of 
approximately n e « 10 11 - 10 13 cm" 3 ). 

In summary, then, standard accretion shock physics does not 
seem to be compatible with a shock-heating scenario for the soft 
spectral component forming the bulk of the Ovn (and Neix) 
lines in T Tau. We emphasize that T Tau is now the second ac- 
creting pre-main sequence star that shows a hig h f/i flux ratio, 
after si milar reports on the Herbig star AB Aur C Telleschi et all 
2006c). This is to be compared with three CTTS previously re- 
ported to show high densities, n e > 10 11 cm" 3 , namely TW Hya 
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dKastner et all 120021), BP Tau dSchmittetalll2005l) . and V4046 
Sgr dGiinther et al.Ll2006l) . 

Abundance anomalies. High abundances of the volatile el- 
ements N and Ne in TW Hya with respect to the Fe abundance 
have been suggested to be due to stronger condensation of Fe 
and some other elements to dust grains in the circumstellar disk 
( Stelz er & Schmittll2004l) . If the star preferentially accretes from 
the gas phase, then an elevated Ne/Fe ratio would be measured 
in the coronal or shock-heated plasma. 

Anomalously high Ne abundances have been found in the X- 
ray source of TW Hya dKastner et all 120021 IStelzer & SchmitiL 
2004), but also in non-accreting evolved stars such as mem- 
bers of the RS CVn class (Ne/Fe abundance ratios up to 13.4, 
Auda rd et alj [2003), and in pre-main sequence stars thought 
to be non-ac creting (weak-line) T Tau stars, such as TWA 5 
(Ne/Fe « 10; [Argiroffi et al. 20051) and HP 98 800 (Ne/Fe * 5; 
iKastner et all r2004bl) . iTelleschi etail d2006bl) have systemati- 
cally analyzed RGS spectra of pre-main sequence stars and com- 
pared Fe/Ne and Fe/O abundance ratios for them. They found 
that i) both ratios tend to be higher for young G-type stars 
than for K- or M-type stars; this trend still holds for near-zero- 
age main-sequence stars; and ii) there is no difference between 
CTTS and WTTS with regard to this trend. 

Fig. [7] shows the Fe/Ne ratio derived from high-resolution 
spectroscopy as a function of spec tral type for T Tau sta rs and 
ZAMS stars (for references, see ITelleschi et all 12006 b). The 
trend for an increasing ratio toward earlier spectral types is 
clearly seen, independently of the evolutionary status (CTTS, 
WTTS, ZAMS). The ratio for T Tau fits in well fo r a K0 star. The 
same holds true for the Fe/O ratio (not shown; see Tellesc hT et alj 
2006b). We conclude that the abundances observed in T Tau fol- 
low a trend that holds for a larger sample of young stars, and this 
trend - as judged from a still small sample of stars - does not 
seem to be related to accretion properties. 

Soft excess. The most outstanding property o f the old CTTS 
TW Hya is its extremely soft X -ray emission dKastner et all 
120021; IStelzer & Schmitli |2004|) . However, no such X-ray 
sources have been iden tified in larg e samples of CTTS in Orion 
dPreibisch et alll2005l) or in Taurus dGudel et alll2006al) . On the 
other hand, ITelleschi et all d2006bl) have suggested that, apart 
from a dominant hot plasma, there is a soft excess in accret- 
ing pre-main sequence stars, best revealed by an unusually large 
ratio between the (summed) photon fluxes in the O vn He-like 
triplet and the O vm Lya line, found to be around unity for CTTS 
but typically less than 0.5 for WTTS. The intrinsic ratio between 
the luminosities of these lines is of course altered by photoelec- 
tric absorption, i.e., the O vn/O vm ratio tends to decrease for in- 
creasi ng A^h- The A^h value f or T Tau compares with the largest 
in the Tellesch i et alj d2006b) sample, yet the observed flux ratio 
is also among the largest. As reported above, we found a ratio of 
1.06 + 0.29. In the unabsorbed model spectrum, the Ovu r line 
is even a factor of two stronger than the O vm Lya line, and in 
fact is the strongest line in the X-ray spectrum. 

T Tau thus confirms the conjecture that CTTS reveal a soft 
excess, despite the dominance of an extremely hard component 
in the overall X-ray spectrum. T Tau is the most extreme case 
for a soft excess in any T Tau star studied by high-resolution 
spectroscopy, except for TW Hya in which the soft component 
entirely dominates. 

Optical and ultraviolet variability. Accretion shocks in- 
duce excess emission at short optical and ultraviolet wave- 
lengths. Variability in these bands seen in accreting T Tau 
stars has therefore been attributed to short-term variations in 
the mass accretion, or simply inhomogeneities at certain Io- 
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Fig. 7. The Fe/Ne abundance ratio in the X-ray sources of 
various pre-main sequence stars and ZAMS stars as a func- 
tion of spectral class. The abundance measurements for stars 
other than T Tau we r e taken from the previous literature and 
from Tellesch i et all d2006bl) (see that paper for a compre- 
hensive description and references; values were converted, if 
necessary, to refer to th e solar pho tospheric abundances o f 
lAnders & Grevessdl 19891 and, for Fe. lGrevesse & Sauvalll 19991) . 
The value determined in the present work for T Tau is shown in 
bold face. Symbols mark different types of stars: squares: ZAMS 
stars; triangles: WTTS; diamonds: CTTS. 



cations in the streams, rather than to flare-like energy release 
events (Gullbring e t all Il996l) . Variability amplitudes are typi- 
cally modest (a few tenths of a magnitude) and last for 0.6 hr 
to several hours. The gas heats onl y moderately, again in con- 
trast to chromospheric gas in flares dGullbring et all ll 996). The 
induced variability may be explained by accretion rate changes 
of the order of 10%, which will then require a scale length of 
the inhomogeneities of 5 x 10 10 cm for the CTTS BP Tau and, 
combined with the fre e fall velocity, correctl y explains the vari- 
ability time scales (see lGullbring et alll 1 9961 for further details). 
It is possible that these inhomogeneities reflect in variations of 
the softest X-ray flux if X-rays are produced in the accretion 
shock, but a dedicated monitoring study of BP Tau with ROSAT 
and ground-based UBVRI photometry revealed an absence of 
correlated events; X-ray variability may thus n ot be related to 
accre tion-induced optical and UV variability (Gullbring et all 
1997). New evidence was provided by the comprehensive XEST 
survey: a comparison of the X-ray light curves of many T Tau 
stars with simultaneous light curves obtained in the U band or a 
shorter-wavelength ultraviolet band showed little correlated be- 
havior, exce pt for a few example s of flares that were observed in 
both bands dAudardet alll2006l) . 

There is mixed evidence in our data. First, both the XMM- 
Newton ultraviolet data and the ground-based U-band data show 
clear evidence for both day-to-day variations and short-term 
variability. The former shows time scales of at least half a day, 
compatible with modulation of accretion hot spots due to rota- 
tion. The short-term events (e.g., at 6 x 10 4 s in Fig. |2l or at 
37.35 d i n the bottom panel of F ig. [3]) are similar to those re- 
ported by Gullbrin g et all 0996). Conclusions are ambiguous, 
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however: During the XMM-Newton observations, there is a long- 
term trend in X-ray hardness that appears to correlate with the 
UV flux; in particular, there is a decay from the beginning of the 
observation that may be related to flare emission. The harder X- 
rays seen in our spectra cannot be produced by accretion shocks 
as the infall velocity is too low (see above). Some of the varia- 
tions of the UV flux may therefore be due to an excess in chro- 
mospheric and transition region heating (producing excess UV 
emission) in concert with coronal heating (producing harder X- 
rays). Also, the short event at 6 x 10 4 s in the UV appears to be 
reflected in the X-ray hardness curve. The temporal relation is 
typical for flares: The ultraviolet excess, a signature of impul- 
sive heating of the chromosphere, precedes the X-ray signal, a 
consequence of filling of magnetic loops with hot plasma. Thus, 
overall a significant contribution to UV variability may come 
from coronal energy release events, although a clear distinction 
is dif ficult. This view is supported by the large survey presented 
bv lAudardet"aT]j2006l) that shows several obvious X-ray flares 
with correlated U-band or ultraviolet features. 

Modulation effects due to rotation are also possible. The 
UBVR observations in Fig.[3]show some evidence (e.g., within 
the middle or the lower panels, or between the upper and middle 
panels) for a correlation between B or V magnitude and the U-B 
color, indicating that the source is bluer when brighter and redder 
when fainter. This is expected for both heated regions and cool 
spots. Heated regions could be due to either accretion hot spots 
or magnetically active chromospheric regions. Ho wever, the in- 
clination angle of T Tau is sma ll, i = 19-23 deg ( I Herbst et all 
U997tlEisloffel & Mundtlll998l) . Rotationally modulated surface 
features must therefore be located at latitudes of approximately 
+20 deg arou nd the equator. The stand ard magnetic accretion 
scenario (e.g.. lCalvet & G ullbring 1998) predicts that accretion 
hot spots form at high latitudes. It is therefore more likely that, 
at least during our UBVR observations, the modulation was 
due to active regions (unless the accretion also occurs along 
more complicated magnetic field structures toward equatorial re- 
gions). The indication of (at least) two minima and two maxima 
in Fig. [3] could be due to two preferred active longitudes pro- 
ducing chromospherica lly active regions near the equator (e.g., 
lBerdvuginaetal]|l998l) . 

7.2. Implications for models 

Our discussion thus points to the following conclusions: i) If 
mass accretion streams induce X-rays in shocks near the surface, 
then densities should be much higher than inferred from the ob- 
served density-sensitive line flux ratios. It is possible that some 
of the observed X-ray emission originates in accretion shocks, 
but only a fraction of the accretion stream, subject to very low 
densities, would be responsible, ii) Abundance anomalies do not 
seem to reflect accretion properties, and T Tau is no exception, 
iii) T Tau is, however, the most outstanding example, apart from 
TW Hya, of an accreting T Tau star with a soft excess in the X- 
ray spectrum. Such excess is not found in non-accreting young 
stars. We thus suggest that the X-ray production of T Tau is in- 
fluenced by the accretion process although the X-rays may not 
form in the bulk of the accretion footpoints. 

Shocks could form in outflows and j ets that are related to the 
accretion process (see iGudef et al. 2005 for the example of DG 
Tau A). In that case, low temperatures and low densities would 
be expected. Since T Tau does drive outflows and is surrounded 
by Herbig-Haro objects (see references in Sect. [TJ, this possi- 
bility cannot be excluded, but explicit evidence is lacking. The 
high-resolution Chandra HRC image (Fig.[TJ does not show any 



evidence for emission that might be related to a jet, in c ontrast to 
what has been found for DG Tau A dGiidelet aUl2005h . The en- 
tire X-ray spectrum can be explained by a multi-thermal plasma 
subject to an overall photoelectric absorption component, again 
in contrast to the sp ectra of proposed jet-related X-ray sources 
(iGiidel et allkooebh . 

An alternative location for the interaction between the 
plasma heating process and accretion is the corona or magneto- 
sphere itself. Accreting material is guided along magnetic field 
lines and can increase the density in these regions considerably 
(similar to the estimates presented for the footpoints of the ac- 
cretion flows, see above). Depending on the stream density and 
the strength of the magnetic fields near the accretion streams, 
the magnetic geometry may be considerably shaped by the mass 
stream itself; magnetic field lines be straightened by the flow, 
perhaps suppressing coronal magnetic energy release so that 
coronal heating may have little or no effect on infalling gas. 

If, however, a small fraction of the mass loads onto magnetic 
fields of active coronal regions, it will lower the overall tempera- 
ture and slightly increase the density there; cooling thus also be- 
comes more efficient. A fraction of the coronal volume affected 
by a small amount of the infalling material could thus produce a 
cool coronal component, i.e., an observa ble soft excess. A sug- 
gestion along those lines was made by iPreibisch et alj |2005) 
who proposed that the lower Lx/^boi ratio of CTTS compared to 
WTTS in Orion is due to accretion flows completely suppressing 
coronal heating in the respective magnetospheric regions. 

We estimate what fraction of the accreting material would 
have to be heated to T — 1.8 MK as observed in the softest com- 
ponent, assuming that the entire soft component is emitted by 
such material (rather than by plasma evaporating from the stel- 
lar surface). We provide this estimate for the limiting case that 
any parcel of inflowing material that becomes subject to coronal 
heating is heated to T at once, after which it is permanently cool- 
ing by radiation, without further heating. The energy rate to heat 
accreting, ionized gas to T — 1.8 MK is E « 3kTN e for hydro- 
gen gas, where N e is the rate of new, cool electrons flowing into 
the heating region; there is an equal rate of protons. Adopting 
Ma CC « (3 - 6) X 10 _8 M o yr~' for the total accretion rate, we find 
N e = (1.1 -2. 2)xl0 42 z s" 1 , wherez is the fraction of the infalling 
gas stream that is heated. Therefore E « (8 - 16) x 10 32 z erg s" 1 . 
We assume that all energy is radiatively lost by the plasma at 
1.8 MK. Modeling the energy losses of the cool plasma compo- 
nent in XSPEC, we find that only 17% of the radiative power is 
emitted above 0.3 keV, or 28% above 0.1 keV. Because the cool 
component emits 5.2 x 10 30 erg s _1 into the 0.3-10 keV band 
(Sect. |6]l, the total radiative energy loss rate of this plasma is 
3 x 10 31 erg s _1 . We conclude that z — 1.9-3.8% of the accretion 
flow would need to be heated to low coronal temperatures. This 
fraction will be smaller if the gas is continually heated while 
radiating. 

We conclude, then, that T Tau implies an important role of 
accretion in the production of high-energy emission. We suggest 
that part of the accreting material adds a cool coronal component 
by being heated to modest temperatures while flowing into ac- 
tive regions. We are not able to predict why the X-ray excess is 
found at temperatures of about 2 MK, but we show that the ad- 
ditional loading of magnetic fields by cool, accreting mass will 
suppress heating to the otherwise common coronal temperatures 
of>10MK. 
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